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Abstract 

The flavor ud and us pseudoscalar correlators are investigated using fami- 
lies of finite energy sum rules (FESR's) known to be very accurately satisfied 
in the isovector vector channel. It is shown that the combination of con- 
straints provided by the full set of these sum rules is sufficiently strong to 
allow determination of both the light quark mass combinations m u + m^, 
m s + m u and the decay constants of the first excited pseudoscalar mesons in 
these channels. The resulting masses and decay constants are also shown to 
produce well-satisfied Borel transformed sum rules, thus providing non-trivial 
constraints on the treatment of direct instanton effects in the FESR analy- 
sis. The values of m u + ma and m s + m u obtained are in good agreement 
with the values implied by recent hadronic r decay analyses and the ratios 
obtained from ChPT. New light quark mass bounds based on FESR's involv- 
ing weight functions which strongly suppress spectral contributions from the 
excited resonance region are also presented. 
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I. INTRODUCTION 



The divergence of the flavor ij axial vector current in QCD is related to the corresponding 
pseudoscalar density by the Ward identity 



(mi + rrij) qii^Qj ■ 



(1) 



As has been long recognized, this fact, together with the analyticity of the correlator, Ilj^g 2 ), 
defined by 



d'xe^ (0|T (0^.(20^/(0)) |0> = (m* + rrijf rL;,(g 2 



(2) 



allows one to write down sum rules which relate the light quark mass combinations mi+rrij to 
the decay constants of the flavor ij pseudoscalar mesons 0. These sum rules, which include 
the basic unsubtracted dispersion relation (involving U"p and/or its derivatives) [0-0], the 
Borel transformed version of this relation |3|,|5|,[7|-[14| , and finite energy sum rules [lC,15-20 



have been used to either place bounds on m u + and m s + m u , or estimate their values. 

The basic forms of these relations are, for the unsubtracted dispersion relation (DR), the 
corresponding Borel sum rule (BSR) M, and finite energy sum rules (FESR's), 
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(4) 
(5) 



respectively, with the spectral function of 11^, s in line 2 of Eq. (|j) the "continuum 
threshold" (beyond which p^ is approximated by its OPE form), M the Borel mass, and 



w 



s) in Eq. (||) any function analytic in the region of the contour. B H'-a (M 2 ) in Eq. 



is the Borel transform of the OPE representation of II^(Q 2 ) |/J. 

The LHS of either Eq. (^) or (|j) can be evaluated using the OPE provided the relevant 
scale (Q or M) is large compared to the QCD scale. For the FESR case, the condition that Sq 
be similarly large is necessary, but not sufficient, to allow reliable evaluation of the LHS using 
the OPE. The reason is that, except at extremely large sq, the OPE is expected to break 
down over some portion of the circle, \s\ = sq, sufficiently near the timelike real axis ||21|| . In 
the flavor ud vector channel, where the spectral function has been determined very accurately 



from hadronic r decay data |22j23l, one can, in fact, verify this breakdown: FESR's involving 



the weights w(s) = s with k = 0,1, 2, 3, which do not suppress contributions from the region 
near the timelike real axis, are typically rather poorly satisfied at scales 2 GeV 2 < Sq < m 2 



At these scales, however, this breakdown turns out to be very closely localized to the 
vicinity of the timelike axis: as soon as one restricts one's attention to weights with even 
a single zero at s — sq, the corresponding FESR's are very accurately satisfied over this 
whole range of sq pi |. Thus, for the "intermediate" scales 2 GeV 2 < sq < 4 GeV 2 which 
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will be of interest to us, we must also include, as a condition for the reliability of the OPE 
representation of the LHS of Eq. (g), the further requirement that w(sq) = 0. We will refer 
to FESR's satisfying this criterion as "pinch-weighted" FESR's (pFESR's) in what follows. 

In the region below s ~ 4 GeV 2 , where the resonances in the channels of interest (ij = 
ud, us) are well-separated, the spectral function will be dominated by contributions from 
the flavor ij pseudoscalar resonances, P. In the convention where f n = 92.4 MeV and 



fx = 113.0 MeV |25|| , the corresponding contribution to pij, ignoring interference, is 

[ Pij (s)] p = 2f 2 P m P B(s) (6) 

where B(s) = 5(s) in the narrow width approximation, with the standard Breit-Wigner 
generalization to finite width, 

B( \ = 1 Tpmp 
{ } 7t[{s-m p y + r p m P ] ' 1 ' 

Experimentally, both f n and fx are very accurately known, while the higher resonance 
(7r(1300) and 7r(1800) for ij = ud and i^(1460) and i^(1830) for ij = us) decay constants 
are unknown at present.^] The positivity of Pij(s), together with the fact that the weights 
appearing in the spectral integrals of Eqs. (§) and @ are > 0, implies that the 7r (or K) 
pole contributions provide lower bounds to these integrals. The same is true for Eq. (Q) as 
long as the weight w(s) employed is positive for < s < s . These lower bounds allow one 
to obtain corresponding lower bounds for m u + and m s + m u [T|] . To actually determine 
m u + rrid and m s + m u , rather than just set bounds on them, however, one must at present 
provide theoretical input for the higher resonance contributions. These contributions cannot 
be expected to be negligible since the fpfnp factors for all P are formally of the same order 
in the chiral expansion. In fact, in existing analyses, the higher resonance contributions 
are typically larger than the ir (or K) pole contributions - as an example, the 7r(1300) and 
7r(1800) contributions to the s°-weighted FESR used to determine m u + m d in Refs. |10|J20[| 
are a factor of ~ 2 — 3 times the ir pole contribution. 

Two approaches to constraining the higher resonance contributions exist in the litera- 
ture. In the first, additional sum rules have been used to provide an estimate of the decay 
constant of the first excited resonance P,^JD] Jl^JT5| -[T8l . In the second, resonance dominance 



has been assumed to be a good approximation, even in the 3n (or Kmr) threshold region, 
and known ChPT expressions for the threshold values of the spectral functions used to 
normalize sums-of-Breit-Wigner ansatze for the higher resonance contributions. Since the 
thresholds are typically several resonance widths (or more) removed from the resonance 
masses, the peak normalizations (the features of the resonance contributions to which the 
sum rule determinations of the m 8 + rrij are dominantly sensitive) will be ambiguous in this 
approach, depending, for example, on the treatment of the s-dependence of the "off-shell 



/tt(1300) an< l /k(1460) could, in principle, be determined using data from hadronic r decay, but 
this would require disentangling these contributions from spin 1 resonance contributions in the 
same region. Neither the ud nor us spin decomposition for the excited resonance region has been 
performed to date. 
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width" . Potential dangers of this threshold normalization approach have been discussed in 
Refs. [^f],^7|]. The situation in the us scalar channel, where the near-threshold behavior 
of the spectral function is significantly constrained by known Kit 1 = 1/2 s— wave phase 
shifts, is particularly instructive. As shown in Ref. [p7[| , the near-threshold spectral function 
implied, through unitarity, by the Ktt phases and the resulting Omnes representation of the 
timelike scalar Ktt cannot be well represented by the tail of a Breit-Wigner resonance form; 
a significant background component, interfering constructively with the resonance contribu- 
tion in the threshold region, is required. The near-threshold normalization of the resonance 
contribution is, therefore, significantly reduced, producing a corresponding reduction in the 
value of the spectral function at the Xq(1430) resonance peak. This reduction is very sig- 
nificant numerically: the Kq(1430) peak value of the us scalar spectral function obtained 
in Ref. p7| (albeit with some additional assumptions about the high-s behavior of the Ktt 
phase and the form of the Omnes representation) is a factor of ~ 3 smaller than that ob- 
tained, using the threshold-resonance-dominance assumption (TRDA), in Ref. [f[T|| . Even if 
one questions the additional assumptions which go into the precise numerical value of the 
reduction in this case, one should bear in mind that the TRDA ansatz for the us scalar 
channel was shown to correspond to a value of the slope of the timelike Kit form factor at 
threshold incompatible with that known from ChPT |]27||. Further evidence of the potential 



problems of the TRDA approach are provided by the results of Ref. [24]. In Ref. [23], the 



TRDA ansatze of Refs. [|10| , P0| for the ud pseudoscalar channel and of Refs. |Tl],0] for the 
us scalar channel were tested using families of pFESR's in which the OPE scales used were 
the same as those employed in the earlier analyses. If the TRDA spectral ansatz for a given 
channel is a good representation of the physical spectral function in that channel, and if 
the scale of the original analysis was such that the OPE representation could be reliably 
employed, then pFESR's constructed using the same spectral ansatz for the same correlator 
should also be well satisfied. It turns out that, in both the ud pseudoscalar and us scalar 
channels, the TRDA ansatz produces a very poor match between the OPE and spectral 
integral sides of the various pFESR's In contrast, the match corresponding to the us 



scalar spectral function of Ref. [27[ is quite reasonable p4j 



In view of the above observations, we do not employ the TRDA ansatz for the ex- 
cited pseudoscalar contributions, but instead constrain these contributions, in analogy 
to the treatment of the isovector vector and scalar channels in Ref. [2J| Q, by ana- 



lyzing simultaneously two continuous families of pFESR's, corresponding to the weights 
w n(v) = (1 — + Ay) and w^(y) = (1 — y) 2 (l + Ay), where y = s/s . As we will show, 
the set of these constraints is sufficiently strong to allow determination of not only the ex- 
cited resonance decay constants, but also the light quark mass combinations. The input 
required for this analysis is outlined briefly in the next section. Our final results, together 
with a discussion of existing quark mass analyses, are provided in Section III while Section 
IV contains our conclusions. 



2 Note that, in the isovector vector channel, if one ignores the experimental spectral data and 
instead uses the pFESR OPE integrals to fit the decay constants of a spectral ansatz consisting 
of a sum of Breit-Wigner resonance contributions, one obtains a value of the p decay constant in 
agreement with the experimental value to better than the experimental error |29[] . 
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II. INPUT FOR THE PFESR ANALYSIS 



The spectral ansatz for the ud pseudoscalar channel is 

Pud (s) = 2flm A J (s - ml) + 2f*rT$B 1 {s) + 2f 2 2 m A 2 B 2 (s) , (8) 
where m 12 are the PDG2000 |25] masses of the 7r(1300) and 7r(1800), /12 are their (as yet 



undetermined) decay constants, and Bi^(s) are the standard Breit-Wigner forms. We have 
employed PDG2000 values for all resonance widths. The corresponding expression for p us (s) 
is obtained by the replacements it -> K, tt(1300) -> if (1460) and tt(1800) -> if (1830). In 
order that this ansatz provide a good representation of the spectrum over the whole range 
required in the pFESR spectral integrals, so cannot be taken much greater than m 2 ,; if 
it is, an unphysical "gap", with little spectral strength, will be present in the integration 
region. We therefore require Sq to remain less than about [m-i + Y 2 ) 2 — 4 GeV 2 . To create 
a good analysis window in sq without at the same time sacrificing good convergence of the 
integrated D = OPE series, we also take so > 3 GeV 2 . 

The ability to avoid unphysical spectral gaps represents a potential advantage of the 
pFESR framework over its BSR counterpart. For BSR's, the continuum threshold, Sq, 
is usually set by requiring an optimal stability window with respect to the Borel mass, 



M. Taking the ij = us analysis of Ref. |12] as an example, and considering the case, 



Aqcd = 380 MeV, which most closely corresponds to the current experimental determination 
of a s (m 2 ), the stability window is optimized for sq between 6 and 8 GeV 2 The resulting 
spectral ansatz, therefore, has a gap with very little spectral strength from about 5 to 
6 or 8 GeV 2 . It is also worth noting that, after Borel transformation, the scale relevant 
to the running coupling in the OPE is ji — M. For the correlators of interest to us the 
convergence of the transformed D = series becomes good only for M 2 greater than about 
2 GeV 2 . Even if one is willing to tolerate a spectral gap by allowing sq ~ 6 GeV 2 , this 
means that Sq/M 2 will be ~ 1 — 2 over much of any putative stability window in M. Such 
a condition signals non-trivial contributions from the "continuum" region, where only a 
relatively crude approximation to the spectral function is being employed. This leaves only 
a small range of M having both good OPE convergence and acceptably small continuum 
contributions (say less than ~ 30% of the D = OPE term). With such a small range of 
M, the BSR constraints are not sufficiently strong to allow a simultaneous determination of 
the quark masses and excited resonance decay constants. In the case of pFESR's, empirical 
evidence from the isovector vector channel suggests that contributions analogous to the 
less reliable continuum BSR contributions (i.e., those contributions from the region of the 
contour \s\ = sq near the timelike real axis, where the OPE is expected to break down) are 
strongly suppressed by the restriction to weights satisfying w(sq) = 0. 

On the OPE side of the sum rules, one must bear in mind that, in scalar and pseudoscalar 
channels, potentially important contributions from direct instantons exist which are not 
incorporated in the standard OPE representation of Ely ]30|]. Such contributions are, in 
fact, needed to produce a Borel transform, B 



n 



ud (M 2 ), which behaves correctly (i.e., is 
independent of M) in the chiral limit [3^-33]. The instanton liquid model (ILM) |34| provides 
a tractable framework for estimating such contributions. In the ILM, an average density 
(related to the value of the gluon condensate) and fixed average size are employed for the 
instanton distribution. Phenomenological constraints require the average instanton size, pi 
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to be ~ 1/0.6 GeV [pT|p2| , p4[ . Instanton contributions to B \n. U d] (M 2 ) then exceed one- 
loop perturbative contributions below M 2 ~ 1 GeV 2 , but drop to less than ~ 15% of this 



contribution for M 2 ~ 2 GeV [ 32| . Q Direct instanton contributions have been neglected in 
recent treatments of the ud and us pseudoscalar channels, apart from the BSR ud analyses 



of Refs. [ pl3||L4|| , both of which employed the ILM. The numerical impact of the neglect of 
these contributions should be small for BSR analyses at scales M 2 > 2 GeV 2 since the Borel 
transform is known to rather strongly suppress ILM contributions with increasing scale?]. 
This is, however, not true of FESR analyses, for which ILM contributions fall off, relative 
to the D = perturbative contributions, much more slowly with increasing so- 

In what follows, we will use the ILM to estimate direct instanton contributions to the 
and pFESR's. ILM contributions to pFESR's corresponding to polynomial weights 



w 



N 



can be evaluated using the result |36 

ds s 
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2m J\ 
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ILM 4ti- 



'° dss^jJp^YAp^) , (9) 



where r] u d = 1, 7] us is an S , f/(3)-breaking factor whose value in the ILM is ~ 0.6 [^2]| , and 
the result is relevant to scales ~ 1 GeV 2 . 

One should bear in mind that phenomenological support for the ILM exists primarily 
for those scales (~ 1 GeV ) where instanton contributions are numerically important in 
pseudoscalar BSR's, and that this scale is significantly lower than that (~ 3 — 4 GeV 2 ) 
relevant to our pFESR analysis. It is, therefore, useful to have an independent test of our 
use of the ILM representation of instanton effects. In this regard, one can take advantage 
of the much stronger suppression of ILM contributions in the BSR framework. The basic 
idea is as follows. One first determines the excited meson decay constants for the channel 
of interest, using the pFESR framework. These values then determine the s < so part 
of the spectral ansatz for a BSR treatment of the same channel. (The spectral function 
for s > Sq is, as usual, approximated using the continuum ansatz; we fix the continuum 
threshold, So, following standard practice, by optimizing the stability of the output, in this 
case, the quark mass combination, rrii + rrij, with respect to M 2 .) For M 2 ~ 2 GeV 2 , 
where (1) convergence of the Borel transformed D = series is still reasonable and (2) 
continuum contributions are still relatively small (not yet exceeding ~ 30% of perturbative 
contributions), the resulting BSR should then allow determination of the only remaining 
unknown, rrii + rrij, with good accuracy. The ILM contributions play little role on the 
OPE+ILM side of the BSR's at these scales, but are important for the pFESR's, and hence 
for the values of the resonance decay constants used as input to the BSR's. If the ILM 
representation of direct instanton effects is reasonable at the scale of the pFESR analysis, 
the pFESR and BSR determinations of rrii+mj, which will then have been obtained using the 



3 The combination of 2-, 3- and 4-loop contributions roughly doubles the Borel transformed 1-loop 
D = contribution at M 2 = 2 GeV 2 , hence further suppressing the ratio of ILM to perturbative 
contributions. 



4 For example, the bound obtained in Ref. [ 14 1 is raised by < 5% if ILM contributions are turned 
off |3l 
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same excited resonance decay constants, should be compatible within their mutual errors. 
Since the continuum approximation for the spectral function is a relatively crude one, and 
the stability criterion for choosing sq typically leaves a gap in the BSR spectral model, there 
are uncertainties in the BSR analysis beyond those associated with the uncertainties in the 
OPE input, which are shared by the pFESR and BSR analyses. In order to get a rough 
estimate of these additional uncertainties we allow so to vary in an interval of size 1 GeV 2 , 
i.e., by ±0.5 GeV 2 about the value corresponding to optimal stability, and assign a ±20% 
error to the size of continuum spectral contributions. Since the Sq values we obtain are 
> 3.7 GeV 2 , we consider the latter estimate sufficiently conservative^. The uncertainties on 
m>i ± rrij induced by use of the continuum approximation are then not large, particularly in 
the region near M 2 = 2 GeV 2 , where the BSR continuum contributions are less than ~ 30% 
of the D = OPE term. The BSR/pFESR cross-check is, as a result, most reliable at these 
scales.^ 

The OPE representation of U'L(Q 2 ) is known up to dimension D = 6, with the dominant 
D = perturbative contribution known to 4-loop order The D = term is given by 

2§ 



n"(Q 2 



3 (mi + rrij 



D=0 



8tt 2 



Q 2 



l 



n 



a + 14.1793a 2 ± 77.3683a 3 



(10) 



where a = a{Q 2 ) = a s (Q 2 )/n, fhk = rrik(Q 2 ), with a s (Q 2 ) and m{Q 2 ) the running coupling 
and running mass at scale fi 2 = Q 2 in the MS scheme. The D = 2 term involves quark 
mass corrections to the leading D = result. For ij = ud it is numerically negligible, while 
for ij = us it is given by [^J 



Ks(Q 2 



D=2 



4vr 2 



m, 



fm 2 



, 28 
l + y a + 



8557 _ 77 
72 3 w 



(11) 



In writing Eq. (|TT|) , we have dropped terms involving m Uj d, except in the overall prefactor 
(rh s + rh u ) 2 . The D = 4 ud contributions are 111 



KAQ 2 



D=4 



{m u + fh d ) (I 4 

\fi, + -an, - 



, 26 

i + T » 



2m < uu > 



28vr 



(12) 



5 For the analogous cases of the ud vector and axial vector channels, where the hadronic spectral 
functions are known experimentally from hadronic r decay data, the maximum deviation of the 
actual spectral function from its 4-loop OPE continuum approximation is less than ~ 1/3 of the 
OPE version in the interval 2 GeV 2 < s < m 2 [22]. Note that these scales are smaller than those 



for which we will be employing the continuum approximation, and that we are concerned with the 
average, rather than maximum, deviation in the range s > sq. 

6 The ratio of the continuum to the D = OPE contribution grows relatively rapidly with M 2 . 
For the ud case, for example, it has already reached ~ 50% by M 2 = 3 GeV 2 and ~ 65% by 
M 2 = 4 GeV 2 . 
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where and flf are the RG invariant modifications of (aG 2 ) and (m s ss) defined in Ref. [11], 
rh = (m u + mrf)/2, and we have dropped numerically negligible terms of 0(m 4 ); the D = 4 
us contributions are, similarly, (Tl[ 



Ks(Q 2 



D=A 




2 < mMU > 



, 23 

1 H a 

3 



(13) 



where we have again dropped terms suppressed by powers of m/m s relative to those shown, 
except in the overall (m s + m u ) 2 prefactor. Finally, the D = 6 contributions are [ I I 



n^(g 2 ) 



(ra s + m u y 



D=e 



32 
~9 



-3 [(rriigqjcr ■ Gqj + m^o ■ Gqi)} 
n-apvsA liqiqi} 2 + {Wj) 2 - 9(qiqi)(q~jqj) 



(14) 



where Pvsa describes the deviation of the four-quark condensates from their vacuum satu- 
ration values. 

Numerical values of the input required on the OPE+ILM side of the sum rules are 
as follows: p r = 1/(0.6 GeV) US], oc s {m 2 T ) = 0.334 ± .022 US ( a ^ 2 ) = (°- 07 ± 
0.01) GeV 4 0, {m u + m d )(uu) = -f 2 ml (the GMOR relation) [] , 0.7 < (ss)/(uu) = 
r c < 1 0,0; (gqaFq) = (0.8 ± 0.2 GeV 2 ) (qq) |9[ and p VSA = 5 ± 5 (z.e., allowing, 
to be conservative, up to an order of magnitude deviation from vacuum saturation for the 
four-quark condensates). The D = 0, 2 and 4 contributions to the OPE integral have been 
evaluated using contour-improvement [^Jill, which is known to improve convergence and 
reduce residual scale dependence pif . For this purpose, we employ the analytic solutions for 
the running coupling and running mass obtained using the known 4-loop-truncated versions 



of the (3 [mi and 7 |43[] functions, with the value of 



noted above as input. 



III. RESULTS AND DISCUSSION 



A. Quark Mass Bounds 



Bounds for the light quark masses based on the known values of the ir or K pole con- 
tributions and the positivity of the spectral function, whether obtained using the dispersion 
formulation, BSR's or FESR's, all depend on the scale employed in the OPE. Since, at the 
scales for which the resulting bounds are of phenomenological interest, the 0(a 2 ) and 0(a 3 ) 
terms in the integrated D = OPE series are not numerically negligible, earlier versions 
of these bounds, based on two-loop and three- loop forms of the D = part of H u d,us, are 



7 Deviations from the GMOR relation have recently been shown to be at most 
resulting error on the m s analysis is completely negligible. 



381. The 



S 



superceded by the work of Ref. || (LRT), which employed the 4-loop OPE expression. The 
bounds of LRT are based on the dispersion relation for H'-p and the higher derivative mo- 
ments thereof. Restricting our attention to the results in LRT corresponding most closely 
to the experimental value a s (m 2 ) = .334, i.e., Aqq D = 380 MeV, the most stringent bounds 
arise from what in LRT is called the "quadratic inequality" ||. These bounds decrease with 
increasing OPE scale, Q 2 , and, for Q 2 = 4 GeV 2 , yield (from Figs. 2 and 3 of LRT) 

[m s + m u ](fi = 2 GeV) > 105 MeV 

[m u + m d ]{ii = 2 GeV) > 8.1 MeV . (15) 

Normally one would expect the convergence of the 4-loop D = OPE series to be quite 
good at scales as large as Q 2 = 4 GeV 2 . In this case, however, the denominator appearing 
on the RHS of the quadratic bound (see Eq. (19) of LRT), which has the form 

-^QCD^QCD _ 2 (fQCD^ =1+ M_ + 61 ?9 _ 2 + 15 _ 3 + . ^ ^ 

is very slowly converging, behaving as 1 + .83 + .61 + .51 at Q 2 = 4 GeV 2 . The bounds 
in Eq. (|T5| ) are thus likely to have a significant residual uncertainty associated with the 
truncation at 0(a 3 ) |. The behavior of the D = series is in fact much better for the 
zeroth moment LRT bound. At the lowest scale shown in Fig. 1 of LRT (Q = 1.4 GeV), 
the behavior of the D = series is 1 + .45 + .22 + .15, already quite well-converged. The 
corresponding bound on m s + m u which, reading from Fig. 1, is 

[m s + mj (fx = 2 GeV) > 80 MeV , (17) 

thus seems to us to be subject to significantly less truncation-induced uncertainty. Although 



the zeroth moment bound for m u + md is not quoted in LRT, the result of Eq. (17), together 



with the result R = 2m s / (m u + mj) = 24.4 ±1.5 determined from ChPT |[44|| , would imply 



[m u + m d ]{n = 2 GeV) > 6.6 MeV . (18) 

The result of Eq. fllTf ) is in good agreement with the bound obtained by the same authors 
from the study of the ud scalar channel^ using constraints on the timelike scalar-isoscalar 
mi form factor from ChPT and tttt phase shift data in the region 4m 2 < s < (500 MeV) 2 
(see Fig. 4 of LRT), 

[m u + m d ](n = 2 GeV) > 6.8 MeV . (19) 



8 If one wished to work, e.g., at a scale such that the 0(a 3 ) term in Eq. ( |l6| ) were less than 
~ 20% of the leading term, one would need to go to Q 2 ~ 9 GeV 2 , at which scale the bounds on 
[m s + m u ](fi = 2 GeV) and [m u + rrid\(Li = 2 GeV) would be reduced to ~ 60 and ~ 3.4 MeV, 
respectively. 

9 The D = OPE series corresponding to this bound has the same (good) convergence behavior 
as that given for the zeroth moment bound above. 



9 



An analogous bound for m s was obtained from a treatment of the us scalar correlator 



employing ChPT constraints for the timelike scalar Kir form factor pq| . Taking the case 
from that reference corresponding to the plausible assumption that the one-loop ChPT 
expression for the Kit form factor is accurate to the .5 — 1% level in the region < s < 
m K ~ m ^ the resulting bound is 

m s (/i = 2 GeV) > 65 MeV , (20) 

which is less stringent than that in Eq. ([T7D . Other recent bounds are (1) that obtained 
in Ref. by combining the upper bound on (qq)(l GeV) allowed by the analysis of the 
D — > K*£u£ vector form factor with the (assumed to be well-satisfied) GMOR relation, 

[m u + m d }(2 GeV) > 6.8 MeV , (21) 



and (2) that obtained in Ref. fllj] using BSR's and Holder inequalities at scales ~ m T 



[m u + m d ](2 GeV) > 4.2 MeV . (22) 

Note that the latter bound was obtained including ILM contributions on the OPE side of 
the sum rule; the bound is ~ 5% higher if ILM contributions are turned off |3"5| . All other 
bounds noted above were obtained neglecting direct instanton contributions. This neglect 
should have little impact on dispersive bounds such as that of Eq. ([17]) since, if one uses the 
ILM to estimate these effects, the lower bound of Eq. ( |TTD is reduced by only 3 MeV. 

An alternate approach to using the positivity of pij to obtain quark mass bounds is to 
employ pFESR's with weights satisfying w(s) > in the region < s < s . A potential 
advantage of this approach is the freedom to choose weights which strongly suppress con- 
tributions from the excited resonance region. Strong suppression of this type should lead to 
bounds which are "close" to the actual mass values. One can arrange such strong suppression 
by choosing w(y) = (1 — y) N p(y) with N sufficiently large. Here p(y) is a "residual polyno- 
mial" which has to be chosen in such a way as to (1) keep the coefficients in w(y) small (thus 
avoiding the growth of unknown higher D contributions) and (2) retain good convergence 
of the integrated D = OPE series. The construction of such weights was considered in a 
different context previously j|7|]. Here we consider quark mass bounds based on pFESR's 



employing the three weights of this type constructed in Ref. ||47||. It turns out that both the 



D = OPE convergence and the stringency of the resulting bounds is best for the case of 
the weight called w 2 o(y) in Ref. |47|], so we present results only for this case. The behavior 
of W2o{y) in the integration region (0 < y < 1) is shown in Fig. 1. (Its explicit form may be 
found in Ref. P7fl -) For so = 4 GeV 2 , the contour-improved D = OPE series for the W20 
pFESR, truncated at 0(a 3 ), converges quite reasonably, behaving as ~ 1 + .55 + .28 + .19. 



10 Without this constraint, working with high powers of the factor (1 —y) typically produces poly- 
nomials with large coefficients for the higher degree y k terms. Since y k terms with large k are 
associated with OPE contributions of large dimension, which are poorly constrained phenomeno- 
logically, large y k coefficients signal potentially large, and essentially unknown, non-perturbative 
contributions 48.49], and hence must be avoided. 
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Moreover, since, for example, if we define yK(iim) = m |;(i46o)/4 GeV 2 , W20 (yK(U60)) = 0.11, 
there will be nearly an order of magnitude suppression of excited resonance contributions, 
relative to the K contribution, in the us channel. Unfortunately, the D = convergence 
deteriorates if one tries to go to lower sq, where this suppression would be much stronger. 
Ignoring possible direct instanton contributions, one obtains 

m s (2 GeV) > 93 MeV . (23) 

The convergence is obviously not sufficiently rapid that one should rule out values of the 
bound a further ~ 5 or so MeV lower. The analogous bound for m u + m d is 

[m u + m d }{2 GeV) > 6.6 MeV . (24) 

These bounds should be compared only to those bounds listed above which also neglect 
possible instanton effects. As expected, the rather strong suppression of excited resonance 
contributions relative to K pole term produces a bound on m s , Eq. (|23|), which is more 
stringent than the zeroth moment LRT bound. The m u + m d bound of Eq. fl2"3|), however, 
remains comparable to the LRT m u + m d bound, though still having the advantage that 
one would expect it to represent a better approximation to the true value. If one now 
incorporates an estimate of direct instanton effects using the ILM, the bounds of Eqs. ( |23|) 
and (|24]) are reduced to 

m s (2 GeV) > 84 MeV (25) 

and 



[m u + m d ){2 GeV) > 5.7 MeV . (26) 

The bound of Eq. fl25|) remains slightly more stringent than that of Eq. ([17]). A more 
stringent bound on m u + m d , 

[m u + m d }(2 GeV) > 6.9 MeV , (27) 



can be obtained using Eq. (|25| ) in combination with the mass ratios obtained from ChPT [31 

To go beyond these bounds, we must attempt to also determine the excited resonance 
decay constants as part of the pFESR analysis. This extension of the analysis is described 
in the next section. 



B. Quark Masses and Excited Meson Decay Constants 



To simultaneously extract rrii + rrij and the corresponding excited pseudoscalar decay 
constants, we perform a combined analysis of pFESR's based on the weight families w^(y) 
and wf)(y) and work in the window 3 GeV 2 < so < 4 GeV 2 . For these so, the D = 
OPE series converges well for all A > 0, and the spectral ansatz should be of the correct 
qualitative form. Larger values of A correspond to larger relative contributions from the 
excited resonance region, and hence are useful for constraining the unknown resonance decay 
constants. To explore sensitivities to the choice of analysis regions, we have also considered 
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the alternate ranges 3.6 GeV 2 < so < 4 GeV 2 and 2 < A < 6, as well as considering 
separate w^(y) and Wp(y) analyses (thus checking the mutual consistency of the pFESR's 
corresponding to the two weight families). The only significant impact of uncertainties in 
the experimental input for the resonance parameters is that occurring in the ud analysis, 
associated with the 7r(1300) width; this is a consequence of the rather wide range, 200 < 
T (7r(1300)) < 600 MeV given in the PDG2000 compilation. In what follows, we quote 
errors from this source separately, labelling them with the subscript £ T". Uncertainties 
associated with changes in the so and A analysis windows and weight family choice are 
added in quadrature and denoted by the subscript "method" . Finally, those errors denoted 
by the subscript "theory" are obtained by combining in quadrature errors associated with 
uncertainties in the OPE input parameters Pvsa, (& S G 2 ), a s (m 2 ) and r c and our estimate 
of the error due to truncation of the dominant D = OPE contribution at 4-loop order. 
The latter is obtained by evaluating the 0(a A ) contribution that would result if we assumed 
continued geometric growth of the coefficients, i.e., the presence of an additional term ~ 
422a 4 in the polynomial factor of Eq. (|10]).[3 It t urns out that, when ILM contributions are 
included, the spectral contribution of the second resonance is small for both channels, and 
hence that the corresponding decay constant can be determined with only limited accuracy. 
When quoting results for the second decay constant in this case, we will, therefore, display 
only the range of values allowed by the combined (i.e., "theory", "method" and (for the ud 
channel only) 'T" ) errors. The analysis of the ud channel has been described briefly already 
in Ref. ||. 

The results obtained from the analysis, when ILM contributions are included on the 
theoretical side of the pFESR's, are as follows. For the ud channel we have 

K + m d ] (2 GeV) = 7.8 ± 0.8 r ± 0.5 theory ± 0A method MeV (28) 
/ ff( i3oo) = 2.20 ± 0.39 r ± 0. 18 theory ± 0. 18 method MeV (29) 
< /.(isoo) < 0.37 MeV , (30) 

and for the us channel 

m s (2 GeV) = 100 ± A theory ± 5 method MeV (31) 
/x(i460) = 21.4 ± l.Qtheory ± 2-3 me thod MeV (32) 
< /^(isao) < 8.9 MeV . (33) 

From Eqs. (|28|) and (^), we see that the uncertainty in the 7r(1300) width is, in fact, the 
dominant source of error in the determination of both m„ + m d and /,r(i3oo) • To get a feel 
for the relative size of the various contributions to the "theory" error we note that, for the 
ud case, the errors in [m u + m d ] (2 GeV) due to the uncertainties noted above on the input 
parameters Pvsa, (&sG 2 ), a s (m 2 ) and truncation at 0(a 3 ) are ±0.25, ±0.05, ±0.28 and 
±0.25 MeV, respectively. The corresponding contributions to the errors on m s (2 GeV) are 
±1.5, ±0.4, ±2.3 and ±3.1 MeV, respectively, with a further contributions of ±0.2 MeV 



In view of the discussion in Section 5 of Ref. [50], this estimate is likely to be a very conservative 
one. 
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due to the range of r c employed in this case. The agreement between the OPE and spectral 
integral sides of the various pFESR's corresponding to the results above is very good. The 
fit quality for the us channel is displayed, for the and families, in Figs. ^| and |3|, 
respectively. The analogous wfj and fits for the ud channel are shown in Figs. 4 and 5 
of Ref. respectively. The ratio R = 25.6 ± 2.6 implied by the above results is in good 



agreement with the value, 24.4 ± 1.5 obtained from ChPT in Ref. | 44| . 

If one repeats the pFESR analysis, but now with the ILM contributions set to zero, one 
finds, for the ud case, 

K + m d ](2 GeV) = 9.9 ± 1.2 r ± 1.0 theory ± 0.5 method MeV (34) 

/tt(1300) — 

2.41 ± 0.50 r ± 0.21 theory ± 0.27 method MeV (35) 

/tt(1800) = 1.36 ± 0.16r ± 0.09^eory ± 0.11 

method 

MeV , (36) 

and, for the us case, 

m s {2 GeV) = 116 ± 1 theory ± ^method MeV (37) 
/if(1460) = 22.9 ± 2.1 t heory ± ^-2 m ethod MeV (38) 
/iqi830) = 14.5 ± l.Btheory ± 0-^method MeV . (39) 

The corresponding OPE/spectral integral match is again excellent. This is illustrated for 
the us case, for the family of pFESR's, in Fig. [|. (The agreement for the corresponding 
u>£ pFESR's as well as that for the ud case is not shown explicitly, but is, in fact, of equal 
quality to that for the us wfj family.) The resulting mass ratio, R = 23.3±2.8, is also in good 
agreement with that obtained from ChPT. We thus see that, while the pFESR fit provides 
a good determination of m 8 + rrij and the resonance decay constants once the form of the 
theoretical side of the sum rule (i.e., whether including or excluding ILM contributions) has 
been fixed, it does not, by itself, provide any additional evidence as to whether inclusion 
or exclusion of these contributions is favored. While inclusion of ILM effects is, of course, 
indicated by arguments external to the pFESR analysis, the pFESR analysis itself shows 
only that, in the absence of these contributions, significantly larger values of the relevant 
quark mass combination and second resonance decay constant are required in both the ud 
and us channels. 

We now turn to the BSR analyses of the ud and us channels, which should provide 
additional constraints on the ILM modelling of instanton effects in the pFESR analyses. 
Expressions for the Borel transforms of the OPE side of the sum rules can be found in 



Refs. |IT| , p!2| , p8| , and that for the Borel transform of the ILM contributions in Ref. [JT3 
We take central values for all OPE input, and employ the corresponding central values for 
the excited resonance decay constants, determined above, as input to the BSR analysis. 
To facilitate the BSR/pFESR comparison, we quote only those errors present in the BSR 
analysis which do not also enter the pFESR analysis, namely those associated with (1) the 
±0.5 GeV 2 variation of the continuum threshold parameter s about its optimal stability 
value and (2) the assumed 20% uncertainty in the size of the continuum spectral contri- 
bution. (Additional errors, associated with uncertainties in the values of the OPE input 
parameters, are common to both analyses, and the corresponding errors, as a result, are 
strongly correlated between the pFESR and BSR treatments.) To be conservative, we take, 
as our estimates for these errors, the maximum change in the value extracted for rrti + rrij 
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in our BSR analysis window (see below) produced by the stated variations in so and the 
magnitude of the continuum contribution. These two sources of error have been combined in 
quadrature in quoting results below. A conventional rule-of-thumb is that the BSR analysis 
window should be restricted to M 2 values for which the perturbative continuum contribution 
is less than ~ 50% of the OPE contribution (for a discussion see, for example, Ref. |p2]| ). 
Since, for the ud case, this corresponds to M 2 less than ~ 3 GeV 2 , we work with a BSR 
analysis window 2 GeV 2 < M 2 < 3 GeV 2 . 

The dependence of [m; + rrij}(2 GeV) on M 2 in the extended range 2 GeV 2 < M 2 < 
4 GeV 2 resulting from the BSR ILM analyses is shown in Fig. |5] for ij = us. (The analogous 



result for ij = ud is shown in Fig. 6 of Ref. |51|.) The solid line corresponds to the 
optimal stability value of so, the upper and lower lines to values 0.5 GeV 2 lower and higher, 
respectively. The quark mass values obtained from this analysis are 

[m u + m d }{2 GeV) = 7.5 ± 0.9 MeV (40) 
m s (2 GeV) = 91 ± 9 MeV . (41) 

These results are to be compared to the central pFESR values of Eqs. ( |2"5| ) and ([TI]) above. 
The consistency of the two determinations is excellent for the ud channel, but only marginally 
acceptable for the us channel. The consistency of the central pFESR and BSR us determi- 
nations can be improved by allowing somewhat larger values of r] us . For example, r] us = 0.8 
produces a central value m s {2 GeV) = 97 MeV, with a corresponding central BSR deter- 
mination 89 ± 9 MeV, while r] us = 1 corresponds to m s (2 GeV) = 92 MeV (pFESR) and 
87 ± 10 MeV. In view of the size of the BSR errors, such improvement cannot be taken as 
physically meaningful; this exercise does, however, indicate that errors comparable in size 
to the difference of the pFESR and BSR central values, associated with the crudeness of the 
ILM representation of instanton effects, may still be present in the pFESR results. We will 
therefore include an additional error, given by this difference of central values, in our final 
version of the errors for the light quark masses. 

For the case that no ILM contributions are included, the BSR results, corresponding 
to central values of all input, and the corresponding central values of the resonance decay 
constants, are 

[m u + m d }{2 GeV) = 8.8 ± 0.6 MeV 

m s (2 GeV) = 100 ± 6 MeV , (42) 

which are to be compared to the central values of Eqs. ([£|) and ([TFp. The pFESR determi- 
nations in both cases lie significantly outside the range allowed by the BSR error. 

Consistency between pFESR and BSR analyses thus favors inclusion of the ILM contri- 
butions. To see that the level of inconsistency between the pFESR and BSR results in the 
absence of ILM contributions is, in fact, significant, the following exercise is useful. Rather 
than optimizing the pFESR analysis by varying simultaneously +rrij, fx and f'2, we may, 
for each value of rrii + rrij, find the values of f\ and f 2 which produce the best OPE/spectral 
integral match. We then use these values of fi and f2, as usual, as input to the correspond- 
ing BSR analysis and look for those values of rrii + rrij for which the pFESR input value is 
compatible with the BSR output value, within the additional errors of the BSR analysis. 

For the ud case, in the absence of ILM contributions, this compatibility is obtained only 
for [m u + m d ]{2 GeV) less than 8.1 MeV (pFESR value)/7.6 MeV (BSR value). Taking the 
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"marginal" case, corresponding to the pFESR value [m u + rrid](2 GeV) = 8.1 MeV to be 
specific, one finds that, although the quality of the OPE+ILM/spectral integral match is 
significantly worse than that for the fully optimized fit above, it is perhaps still acceptable 
(see Fig. || for the fit quality for the wd case; the quality is comparable, though marginally 
better, for the wn case). Thus, in this case, although the inclusion of ILM contributions is 
favored, we do not consider it possible to rule out their absence. Note, however, that the 
analysis, in the absence of ILM contributions, is only self-consistent for values of m u + 
compatible with those obtained from the analysis including ILM contributions. The value 
of /,r(i3oo) obtained in this case, 1.74 MeV, also turns out to be compatible, within errors, 
with that given by Eq. (^). 

For the us case, in the absence of ILM contributions, compatibility is achieved only 
for m s {2 GeV) less than 94 MeV (pFESR value) /89 MeV (BSR value). The "best" fit 
pFESR solution for such a value of m s , however, represents an extremely poor quality 
OPE+ILM/spectral integral matchP]. We thus find no acceptable, consistent spectral solu- 
tion in the us case without the inclusion of ILM contributions. This, of course, also favors 
the inclusion of such contributions for the ud channel. 

In view of these observations, we take as our final results those corresponding to the 
pFESR analysis with direct instanton contributions estimated using the ILM. Including our 
additional estimate of the error associated with the crudeness of the ILM, and combining 
all sources of error in quadrature, our final results for the light quark masses become 

[m u + m d ](2 GeV) = 7.8 ± 1.1 MeV (43) 
m s (2 GeV) = 100 ± 12 MeV . (44) 

Since the pFESR/BSR consistency is excellent for the ud channel, but marginal for the us 
channel, an alternate determination of m s , using the result of Eq. fl43|) above in combination 
with the ChPT-determined mass ratio R = 24.4 ± 1.5, might be preferable. The result of 
this determination, 

m s (2 GeV) = 95 ± 15 MeV , (45) 

is in good agreement with that of Eq. (fPj), with only slightly larger errors. Recall that 
self- consistent versions of the combined pFESR/BSR analysis in which direct instanton 
contributions are neglected, in fact, yield values for the light quark masses completely com- 
patible with those of Eqs. fl43|) and fl44|) . For the resonance decay constants, we note that, 
although the value of the second resonance decay constant is sensitive to whether or not one 
includes ILM contributions, that of the first resonance is largely insensitive to the presence 
or absence of ILM contributions, the central values differing by considerably less than the 
uncertainties on the individual determinations. We thus believe that, although the ILM 
may represent a relatively crude model for implementing direct instanton effects, the deter- 
mination of the 7r(1300) and K(1460) decay constants given by Eqs. ( p9~D and (|32| ) should 



12 The average OPE/spectral integral discrepancy over the so,A analysis window is, for example, 
23% for the w D pFESR family. 
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be reliable to within the stated errors. Combining these errors in quadrature we then have, 
for our final results, 

U {im = 2.20 ± 0.46 MeV (46) 
and /if(H60) = 21.4 ±2.8 MeV . (47) 

That these values differ by a factor of ~ 10 is compatible with the fact that the excited 
pseudoscalar decay constants vanish in the chiral limit, and hence are proportional to the 
relevant quark mass combination near that limit. 



C. Discussion 



Other recent sum rule analyses exist for the pseudoscalar ud |TTTjj2rj] and us |il| , |12f chan- 



nels. In addition, sum rule analyses of the us scalar correlator [11,27,28,53,54], and of flavor 



breaking in hadronic r decay [p^, ^ , p5| -|58| , have been used to extract m s . In this section 
we discuss the relation of our work to that of these earlier references. 

For the ud pseudoscalar channel, Ref. |20[ (P98) represents an update of Ref. pi. (The 
latter employed 3-loop versions for the OPE D = contribution, the running mass and 
the running coupling, P98 the 4- loop versions). We therefore restrict our discussion to the 
latter analysis. The resonance part of the P98 spectral function is of the TRDA form, but 
rescaled by an overall factor 1.5. Two points should be borne in mind regarding the value 
quoted for m u + m d in P98. The first is that the analysis is based on FESR's involving 
the weights w(s) = 1 and s. For these weights, however, the corresponding vector isovector 
channel FESR's are not well-satisfied at the scales employed in P98. (The OPE side has a 
significantly weaker s dependence than the spectral integral side, the latter being obtained, 
in this case, from experimental r decay data ||24||.) The second point is that the ratio of 



quoted values for the running mass at scales 1 and 2 GeV, m(l GeV)/m(2 GeV) = 1.31 ||20 
differs from that, 1.38, obtained using 4-loop running with the central ALEPH determination 
of a s {m T ) as input. The results of P98 thus correspond to a smaller value, a s {m T ) = 0.307, 
the effect of which would be to produce a larger value of m u + m d . One would thus expect 
a poor match between the OPE and spectral integral sides of pFESR's employing the P98 
spectral ansatz and central m u + m d value in combination with current central values for 
the OPE input. This is confirmed by the results shown in Figs. [7] and which correspond, 
respectively, to the output from the wfj and pFESR weight family analyses, in our s , A 
analysis window, obtained using the P98 central value for m u + m d and the P98 spectral 
ansatz. If one performs a re-analysis, still using the P98 spectral ansatz, but now optimizing 
the value of m u + m d using the pFESR approach, one finds, using central values for all OPE 
input, and including ILM contributions, 

[m u + m d ](2 GeV) = 6.8 MeV . (48) 

The same analysis, without ILM contributions, similarly yields 

[m u + m d ](2 GeV) = 7.3 MeV . (49) 

Both of these values are, in fact, in agreement with those corresponding to the upper part 
of the Sq range displayed in Fig. 2 of P98, though not with those for s ~ 2 GeV 2 . In both 
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cases, however, the quality of the OPE(+ILM)/spectral integral match is much inferior to 
that obtained obtained using the solutions for m u + m d , /^(i3oo) an d /^(isoo) above. The 
optimized match is significantly better when ILM contributions are included than when 
they are not. However, in spite of optimization, the consistency between the wfj and 
pFESR families is not good for the P98 spectral ansatz: as shown in ||51||, the match for wfj 



is best where that for is worst, and vice versa (see Figs. 2,3 in |5l| ). 

For the us pseudoscalar channel, the BSR analyses of Refs. JIT] ( JM) and |T2j (DPS) both 
employ a TRDA construction for the if (1460) and if (1830) contributions to the spectral 
ansatz, but differ in their assumptions about the relative sizes of the two resonance decay 
constants: JM assume ' f\m\ = 0.25, DPS that the spectral contributions of the 

two resonances at threshold are approximately equal (for PDG2000 values of the masses 
and widths, this corresponds to f^m^j f\m\ ~ 1.8). The two analyses also differ in their 
treatment of the theoretical side, JM employing 3-loop expressions on the OPE side and 
DPS the 4-loop expressions which become available subsequent to the publication of the JM 
paper. 

We have updated the JM BSR analysis to include 4-loop contributions to the running 
mass, coupling and D = OPE term. For OPE input we use the values employed in our 
analyses above. Including ILM contributions on the theoretical side of the BSR, we then 
find that the JM spectral ansatz, Pjm, corresponds to 

m s {2 GeV) = 96 ± 7 MeV . (50) 

Neglecting instanton contributions, as in JM, we obtain instead 

m s (2 GeV) = 98 ± 6 MeV . (51) 

(The errors in these equations have the same meaning as those for the BSR analyses above.) 
If, however, we employ Pjm, as input, not to a BSR analysis, but to our usual pFESR 
analysis, we find for our central values 

m s (2 GeV) = 107 MeV (52) 

if ILM contributions are included, and 

m s (2 GeV) = 111 MeV (53) 

if they are not. The fit quality for the optimized match is rather poor when ILM contributions 
are included, but is quite good when they are not. The latter point is illustrated for the 
wfj family of pFESR's in Fig. ^| (the quality of the match for the family, which is not 
shown, is even better). Despite the existence of both a good quality pFESR OPE/spectral 
integral match and an excellent BSR stability window, however, we see that the no-ILM 
pFESR and BSR m s determinations based on pju are inconsistent, just as was the case for 
the determinations associated with the spectral ansatz based on the values of fK(U60) and 
(1830) obtained from the no-ILM pFESR analysis. This is, in fact, not surprising, since 
the optimized pFESR spectral ansatz turns out to be rather similar to pjm, the if (1460) 
decay constants of the two models, for example, differing by less than 6%. 

In discussing the DPS analysis of the us pseudoscalar channel, one should bear in mind 
that the result quoted by DPS, m s (l GeV) = 155 ± 25 MeV corresponds to (1) an average 
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of the values obtained using 

A.qcd = 280 MeV and 380 MeV, (2) an average over values 
associated with a range of s , and (3) neglect of m u in the OPE prefactoiQ Since the choice 
A.q CD = 280 MeV is not consistent with the ALEPH determination of a s (m T ), we restrict 

(3) 

our attention to the DPS results obtained using Aq CD = 380 MeV, which corresponds very 
closely to the central ALEPH determination. Restoring m u in the overall OPE prefactor, 
and reading off from Fig. 2 of DPS, concentrating on the curve corresponding to Sq = 
6 GeV 2 , which displays the best stability of m s with respect to M 2 , the central DPS BSR 
determination becomes m s (2 GeV) = 97 MeV. Since the details of the spectral ansatz 
employed are not fully specified in DPS, we are unable to quote errors equivalent to those of 
our BSR analyses above. If, however, we fix the ratio of decay constants in such a way as to 
ensure exact equality of the i^(1460) and fC(1830) contributions to the spectral function at 
physical threshold, and neglect ILM contributions, as in DPS, we find that, after performing 
our usual pFESR analysis, the resulting spectral ansatz, run through a BSR analysis with 
Sq = 6 GeV 2 , reproduces the DPS central value exactly. Estimating our BSR errors as for 
the analyses above, we then have, for our DPS-like BSR determination, 

m s (2 GeV) = 97 ± 6 MeV . (54) 

The pFESR OPE/spectral integral match corresponding to this BSR determination is rea- 
sonable (see Fig. for the wfj family case; the fit quality for the family is not shown, 
but is better than for the family). The central no-ILM pFESR m s value, 

m s (2 GeV) = 109 MeV , (55) 

however, is again inconsistent with the corresponding BSR value. The situation is not 
improved by including ILM contributions: re-doing the pFESR analysis, still with the con- 
strained form of the spectral ansatz, but now incorporating ILM contributions on the theo- 
retical side, one finds a poor quality optimized OPE+ILM/spectral integral match. 

We conclude this section with a reminder of the values obtained for m s via sum rule 
analyses of other channels. Recent treatments of the correlator of <9 M (s'y^u) , for 

which the low-s part of the spectral function is constrained by Kn phases, yield values of 
m s (2 GeV) in the range 115 ± 25 MeV, compatible with either the ILM or no-ILM results 
above. Assumptions about the form of the Omnes representation of the timelike scalar Ktc 
form factor, and the behavior of the Kn phase in the region s > 2.9 GeV 2 , where experi- 
mental phase data does not exist, however, enter the construction of the spectral function 
used in those analyses, so that a significant theoretical systematic error is present, in addi- 



tion to the errors quoted in Refs. p7| , |53| , p4| . A much cleaner approach, in principle, is the 
extraction of m s via pFESR analyses of the flavor-breaking difference of ud and us vector- 
plus-axial-vector correlator sums. The hadronic spectral function required in this case is 
measurable in hadronic r decay. There are two basic complications, first, that the OPE 
representation of the longitudinal contribution to the r hadronic decay width is very badly 



Restoring m u to the prefactor, using ChPT values for the quark mass ratios, and converting to 
the scale u = 2 GeV, the DPS result becomes m s (2 GeV) = 109 ± 18 MeV. 
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behaved at those scales which are kinematically allowed |59|-|62| and, second, that, because 



of the rather strong cancellation in the ud-us spectral difference, the extracted value of 
m s is quite sensitive to even ~ 1% uncertainties in the value of \V US \. The first problem 



can be handled by appropriate weight choices [47|]. The second is numerically relevant be- 
cause the central values of the determinations of \V US \ based on (1) experimental Kg 3 data, 
\V US \ = 0.2196 ± 0.0023 and (2) CKM unitarity, in combination with the experimental value 
of \V u d\, \V US \ = 0.2225 ± 0.0035Q, while consistent within errors, differ by ~ 1.3%. There 
has also been some confusion in the literature resulting from the use, in the various recent 
theoretical analyses, of three different sets of values for the weighted ud-us spectral differ- 
ences, corresponding to three different values of B us , the total (V + A) branching fraction 
into strange hadronic states^ . The strong ud-us cancellation makes the extracted value 
of m s quite sensitive to the (apparently rather small) differences between these B us values. 
The discrepancies between the various values of m s reported in the literature, all of which 
are nominally based on the "same" (ALEPH) r decay data, turn out to be almost entirely 



a reflection of this sensitivity. The situation is discussed in some detail in Ref. ||64|| , where 
the various analyses have also been updated to reflect the current experimental situation (as 
reported in Ref. |[58|| ). Once common input is employed, all hadronic r determinations of 



m s are in excellent agreement |54J . The dominant uncertainty remains that associated with 
\V US \. Using central values of \V u d\ and \V US \ corresponding to either (1) the PDG2000 best 
independent individual determinations (CKMN) (|T4d| = 0.9735 and \V US \ = 0.2196) or (2) 
the PDG2000 unitarity-constrained fit (CKMU) (\V ud \ = 0.9749 and |K S | = 0.2225), one 
obtains 

m s (2 GeV) = 101 ± 18 MeV (CKMN) , (56) 

and 

m s (2 GeV) = 114 ± 16 MeV (CKMU) (57) 
respectively []64| . Either of these results is compatible with that obtained from the pseu- 



doscalar channel analyses above. 



IV. CONCLUSIONS 



We have determined m u + rrid, m s , and the decay constants of the 7r(1300) and 7^(1460) 
with good accuracy from a combined pFESR/BSR study of the ud and us pseudoscalar cor- 
relators. Our results show that it is important to require the consistency of the two different 



14 The central value and "errors" quoted here correspond to the mid-point and extent of the 
PDG2000 unitarity-constrained fit range. 

15 These three values, which are in the ratios 1 : 1.04 : 1.05 correspond to (1) the preliminary (1998) 



ALEPH analysis of strange decay modes |63| |, (2) the final (1999) version of this analysis [55], and 
(3) the recent update (2000) reported by Davier |5S|| . Larger values of B us correspond to smaller 
values of the ud-us difference, and hence to lower values of m s . 
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sum rule approaches. Indeed, we have seen that there exist ansatze for the hadronic spec- 
tral functions which produce both extremely good BSR stability plateaus and high-quality 
pFESR OPE(+ILM)/spectral integral matches, but for which the output quark mass combi- 
nations are inconsistent. This means that BSR or pFESR treatments, by themselves, do not 
provide sufficiently strong constraints to allow one to simultaneously constrain the unknown 
quark masses, unknown resonance decay constants and the theoretical modelling of direct 
instanton effects. The combination of the two approaches does, however, provide sufficiently 
strong constraints. The consistency of the combined analysis is particularly compelling for 
the ud case. The values obtained for the light quark masses are in excellent agreement 
with determinations from other sources, giving us further confidence in the reliability of the 
combined analysis. The corresponding determinations of the 7r(1300) and X(1460) decay 
constants are accurate to 20% and 10% respectively. The latter determination is relevant to 
future improvements in the extraction of m s from hadronic r decay data. While B factory 
data will dramatically reduce the errors on the experimental us vector-plus-axial-vector r 
decay distribution, the ability to use this improvement to significantly reduce the errors on 
the corresponding determination of m s will depend on one's ability to work with pFESR's 
involving weights for which the ud-us cancellation is significantly reduced, so that the errors 
resulting from the uncertainty in \V US \ will, as a result, play a significantly reduced role. 
The existence of significant theoretical systematic uncertainties in versions of this analysis 
which include longitudinal OPE contributions |62| means that "non-inclusive" analyses (in- 
volving only the sum of spin and 1 correlator components) will eventually be required. A 
knowledge of the decay constants of the excited strange pseudoscalar and scalar resonances 
allows a straightforward subtraction of the longitudinal contributions to the experimental 
distributions. In the absence of an experimental spin separation, sum rule determinations 
of the strange scalar and pseudoscalar resonance decay constants with an accuracy even a 
factor of three worse than that obtained above for the A^(1460) are already extremely useful 
as input for such non-inclusive analyses. 
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FIGURES 
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FIG. 1. The behavior of the weight u>2o(y) i n the pFESR integration region. 
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FIG. 2. The OPE+ILM versus hadronic (spectral integral) sides of the us family of 
pFESR's, for m s + m u , fK(i460) an d /k(is30) given by the central values of Eqs. (|3l|), ( |32"| ) and (|33|). 
The solid lines are the hadronic integrals, the dashed lines the corresponding OPE integrals. The 
lower, middle and upper lines in each case correspond to A = 0, 2 and 4, respectively. 
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FIG. 3. The OPE+ILM versus hadronic (spectral integral) sides of the us family of pFESR's 
for m s + m u , fK(l480) an d /ftT(l830) given by the central values of Eqs. (fy), (^) and (|3^). The 
identification of OPE and hadronic integrals, and the cases A = 0, 2, 4 is as for Fig. [2] above. 
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FIG. 4. The OPE versus hadronic (spectral integral) sides of the us family of pFESR's for 
m s + m u , /x(i460) an d /ftr(i830) given by the central values of Eqs. (|37|), ( |38|) and (|^), i.e., in the 
absence of ILM contributions. The identification of OPE and hadronic integrals, and the cases 
A = 0, 2, 4 is as for Fig. above. 
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FIG. 5. The value of [m s + m u ](2 GeV), as a function of the square of the Borel mass, M 2 , 
extracted from the BSR analysis of the us pseudoscalar correlator described in the text. The solid 
line corresponds to so = 4.22 GeV 2 , which produces optimal stability for m s + m u with respect 
to M 2 in the window 2 GeV 2 < M 2 < 3 GeV 2 . The lower (short) dashed line corresponds to 
s = 4.72 GeV 2 and the upper (long) dashed line to s = 3.72 GeV 2 . 
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FIG. 6. The ud OPE/spectral integral match obtained for the pFESR family using as 
pFESR input the value [m u + m d ](2 GeV) = 8.1 MeV, the largest pFESR input for which pFESR 
and BSR values of m u + md are consistent. All notation is as for the pFESR figures above. This 
largest "marginal" m u + value produces the best OPE/spectral integral match among those 
input values for which the pFESR input and BSR output values are consistent; the fit quality, 
moreover, deteriorates rapidly as one goes to lower values of the pFESR input. 
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FIG. 7. The ud OPE/spectral integral match obtained for the pFESR family using the 
central values of all OPE input, the quoted P98 value of [m u + m,d](l GeV) and the P98 spectral 
ansatz. All notation is as for the pFESR figures above. 
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s (GeV 2 ) 

FIG. 8. The ud OPE/spectral integral match obtained for the pFESR family using the 
central values of all OPE input, the quoted P98 value of [m u + md](l GeV) and the P98 spectral 
ansatz. All notation is as for the pFESR figures above. 
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FIG. 9. The us pseudoscalar OPE/spectral integral match obtained for the pFESR family 
using the JM spectral ansatz, central values of all OPE input, and no ILM contributions, after 
optimization of m s + m u in a combined w^, pFESR analysis. All notation is as for the pFESR 
figures above. 
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FIG. 10. The us pseudoscalar OPE/spectral integral match for the pFESR family involving 
the spectral ansatz obtained from a combined w^, pFESR analysis after imposing the DPS-likc 
constraint on the ratio of iT(1460) and -fT(1830) decay constants. The results correspond to central 
values of all OPE input, and to neglect of ILM contributions. All notation is as for the pFESR 
figures above. 
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